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SUMMARY 


The  work  summarized  in  this  report  covers  the  fourth  quarter  of  a  program  with 
a  goal  that  is  twofold:  first,  to  study  the  properties  of  native  point  defect  in  infrared 
focal-plan  array  (IRFPA)  active  and  substrate  materials,  and  second,  to  study  the 
properties  of  native  point  defects  in  two  classes  of  photonic  materials,  the  wide-gap 
II- VI  compounds  (ZnSe  as  the  prototype  for  which  impurity  properties  will  also  be 
calculated)  and  the  nonlinear  optical  materials  (LiNbOa  as  the  prototype).  Our 
accomplishments  in  the  fourth  quarter  include 

•  Completion  of  a  study  of  the  gradient  correction  to  the  local  density 
approximation,  which  is  an  essential  ingredient  to  our  calculations  of 
accurate  absolute  defect  concentrations  in  HgCdTe,  ZnSe,  and  LiNbOa. 

•  Calculation  of  the  degeneracy  factors  for  the  the  single  and  double 
donor  and  acceptor  states  of  the  native  defects. 

•  Analysis  of  strategies  for  calculating  the  localized  states  in  ZnSe. 

•  Development  of  a  preliminary  description  of  the  tellurium  inclusion  for¬ 
mation  process  and  their  native  point  defect  atmospheres  in  HgCdTe, 
and  proposal  of  a  mechanism  for  the  mobility  reduction  in  mercury- 
annealed  liquid  phase  epitaxy  (LPE)  material. 

•  Development  of  input  for  defect  calculations  in  LiNbOa,  and  beginning 
of  calculation  of  the  lithium  vacancy  formation  energy  in  LiNbOa. 

•  Calculation  of  the  electron  mobility  in  Hgo.32Cdo.7sTe  from  a  solution  of 
the  full  Boltzmann  transport  equation  with  Fermi-Dirac  statistics  and 
an  accurate  band  structure.  Excellent  agreement  with  experiments  has 
been  obtained 

•  Calculation  of  transport-related  properties  (Hall  coefficient,  Fermi  level, 
and  energy  gap)  using  accurate  band  structures  and  Fermi-Dirac  statis¬ 
tics.  Results  differ  appreciably  from  those  using  traditional  approxima¬ 
tions. 
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1.  NATIVE  POINT  DEFECTS  IN  HgCdTe.  ZnSe  AND  LiNb03 


1.1  GRADIENT  CORRECTION  TO  THE  LOCAL  DENSITY  FORMALISM 

During  this  quarter,  we  completed  an  important  component  of  our  work  on  defects 
in  HgCdTe,  ZnSe,  and  LiNb03,  which  is  the  gradient  correction1  to  the  local  density 
(LD)  formalism.  The  gradient  correction  (GC)  is  important  in  giving  our  theoretical 
numbers  the  accuracy  necessary  to  calculate  absolute  defect  concentrations,  a  task 
that  to  our  knowledge  no  others  have  attempted  to  date. 

Although  we  have  been  using  the  gradient  correction  for  some  time,  a  recent 
paper  in  the  literature2  has  raised  some  issues  about  the  systematics  of  the  correction. 
We  have  thus  felt  it  necessary  to  undertake  a  more  careful  analysis  of  the  gradient 
correction,  both  to  address  the  issues  raised  by  the  paper’s  authors,  for  our  own 
understanding,  and  to  establish  the  reliability  of  our  defect  work,  which  will  include 
the  correction. 

In  Tables  1  through  4  we  show  the  results  of  our  study  of  the  gradient  correction 
on  the  lattice  constant,  bond  energy,  bulk  modulus,  and  cn— cn  elastic  constant, 
respectively.  Shown  in  these  tables  are  the  results  for  the  diamond  and  zincblende 
semiconductors.  Below  we  discuss  in  turn  the  results  in  each  table. 

As  can  be  seen  from  Table  1,  the  lattice  constants  calculated  with  the  LD  and 
GC  are  both  in  good  agreement  with  the  experimental  results.  The  GC  consistently 
lengthens  the  lattice  constant  with  respect  to  the  LD;  in  some  cases  this  improves 
the  agreement  with  experiment,  in  others  it  worsens  the  agreement.  In  both  cases 
the  average  agreement  with  experiment  of  this  ab  initio  theory  is  better  than  <  1%, 
which  is  quite  remarkable. 

The  bond  energy  is  defined  as  the  magnitude  of  the  crystal  energy,  per  bond, 
relative  to  the  free  atom  energies  of  the  constituents.  For  example  the  bond  energy, 
Eb  for  ZnSe  is  defined  via: 


AEb  +  ZnSetoj,^  -♦  Zn/ree  atom  +  Se/  ree  atom 

The  bond  energy  is  multiplied  by  four  in  this  reaction  because  four  ZnSe  bonds  are 
broken  in  the  reaction.  With  one  exception  (InP),  the  agreement  of  the  bond  energy 
with  the  experimental  values  is  improved  significantly  when  the  gradient  correction 
is  included  in  the  calculation.  In  all  cases  the  LD  overestimates  the  binding  in  the 
semiconductors  -  that  is,  the  cohesive  energies  are  too  large.  The  GC  bond  energies 
are  reduced  with  respect  to  the  LD;  in  some  cases  this  results  in  a  slight  underbinding. 
If  we  consider  HgTe  (the  primary  constituent  of  low  x  Hgi_*Cd*Te  and  ZnSe),  the 
improvement  in  the  bond  energies  is  dramatic:  for  HgTe  we  go  from  a  35%  overbinding 
in  LD  to  only  a  6%  underbinding  with  the  GC.  For  ZnSe,  the  LD  predicts  a  bond 
energy  that  is  27%  too  large,  while  the  GC  predicts  only  a  5%  overbinding.  Because 
we  are  interested  in  calculating  total  energies  in  which  the  free  atom  is  the  reference 
state,  for  example  the  formation  of  a  Zn  vacancy  via  the  reaction 
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Table  1.  Comparison  of  lattice  constant,  in  atomic  units,  for  various  zincblende  semi¬ 
conductors  using  the  local  density  (LD)  approximation  and  the  gradient  correction 
(GC)  to  the  local  density.  Shown  for  comparison  are  the  experimental  values.  Values 
shown  in  parentheses  are  the  percentage  difference  from  the  experimental  values. 


Compound 

Experiment 

LD 

GC 

Si 

10.27 

10.21(-1%) 

10.28(<  1%) 

Ge 

10.69 

10.65(<  1%) 

10.74(<  1%) 

A1P 

10.33 

10.27(— 1%) 

10.36(<  1%) 

AlAs 

10.70 

10.65(<  1%) 

10.73(<  1%) 

AlSb 

11.60 

11.57(<  1%) 

11.68(4-1%) 

GaP 

10.30 

10.21(-1%) 

10.30(<  1%) 

GaAs 

10.69 

10.63(-1%) 

10.73(<  1%) 

GaSb 

11.50 

11.50(<  1%) 

11.62(+1%) 

InP 

11.09 

11.03(— 1%) 

11.13(<  1%) 

InAs 

11.45 

11.43(<  1%) 

11.54(4-1%) 

InSb 

12.25 

12.25(<  1%) 

12.38(4-1%) 

ZnS 

10.23 

10.02( — 2%) 

10.16( — 1%) 

ZnSe 

10.71 

10.55(— 1%) 

10.68(<  1%) 

ZnTe 

11.53 

11.38(-1%) 

11.51(<  1%) 

CdS 

11.00 

10.89( — 1%) 

11.03(<  1%) 

CdSe 

11.44 

11.38(— 1%) 

11.53(4-1%) 

CdTe 

12.25 

12.16(— 1%) 

12.32(4-1%) 

HgS 

11.06 

10.94(— 1%) 

11.09(<  1%) 

HgSe 

11.50 

11.46(<  1%) 

11.60(4-1%) 

HgTe 

12.21 

12.20(<  1%) 

12.33(4-1%) 

Table  2.  Comparison  of  bond  energy,  in  electron  volts,  for  various  zincblende  semi¬ 
conductors  using  the  local  density  (LD)  approximation  and  the  gradient  correction 
(GC)  to  the  local  density.  Shown  for  comparison  are  the  experimental  values.  Values 
shown  in  parentheses  are  the  percentage  difference  from  the  experimental  values. 


Compound 

Experiment 

LD 

GC 

Si 

Ge 

2.58(411%) 

2.22(414%) 

2.31(  <  1%) 
1.88(  -3%) 

A1P 

2.13 

2.35(410%) 

2.05(  -4%) 

AlAs 

1.89 

2.17(4-15%) 

1.85(  -2%) 

AlSb 

1.76 

1.91(  +9%) 

1.61(  -9%) 

GaP 

1.78 

2.08(417%) 

1.76(  -1%) 

GaAs 

1.63 

1.91(417%) 

1.56(  -4%) 

GaSb 

1.48 

1.70(415%) 

1.36(  -8%) 

InP 

1.74 

1.89(  49%) 

1.56(— 10%) 

InAs 

1.55 

1.77(414%) 

1.42(  -8%) 

InSb 

1.40 

1.60(414%) 

1.26(— 10%) 

ZnS 

1.59 

1.82(414%) 

1.53(  -4%) 

ZnSe 

1.29 

1.64(427%) 

1.35(  45%) 

ZnTe 

1.20 

1.43(419%) 

1.15(  -4%) 

CdS 

1.42 

1.63(415%) 

1.33(  -6%) 

CdSe 

1.21 

1.49(423%) 

1.19(  -2%) 

CdTe 

1.10 

1.33(421%) 

1.04(  -5%) 

HgS 

1.02 

1.29(426%) 

0.94(  -8%) 

HgSe 

0.85 

1.19(440%) 

0.84(  -1%) 

HgTe 

0.81 

1.09(435%) 

0.76(  -6%) 

Table  3.  Comparison  of  bulk  modulus,  in  units  of  10nergs/cm3,  for  various  zincblende 
semiconductors  using  the  local  density  (LD)  approximation  and  the  gradient  correc¬ 
tion  (GC)  to  the  local  density.  Shown  for  comparison  are  the  experimental  values. 
Values  shown  in  parentheses  are  the  percentage  difference  from  the  experimental  val¬ 
ues. 


Compound 

Experiment 

LD 

GC 

Si 

9.92 

9.61(  -3%) 

9.32(  -6%) 

Ge 

7.65 

7.13(  -7%) 

6.77(— 11%) 

AlSb 

5.82 

5.58(  -4%) 

5.32(  -9%) 

GaP 

9.14 

8.94(  -2%) 

8.50(  -7%) 

GaAs 

7.84 

7.30(  -7%) 

6.86( — 12%) 

GaSb 

5.79 

5.38(  -7%) 

4.97(— 14%) 

InP 

7.25 

7.01(  -3%) 

6.61(  -9%) 

InAs 

5.79 

5.92(  +2%) 

5.52(  -5%) 

InSb 

4.82 

4.59(  -5%) 

4.22(— 12%) 

ZnS 

7.64 

8.57(+12%) 

7.87(  +3%) 

ZnSe 

6.46 

7.13(+10%) 

6.42(  -1%) 

ZnTe 

5.28 

5.49(  +4%) 

4.98(  -6%) 

CdTe 

4.21 

4.58(  +9%) 

4.05(  -4%) 

HgS 

6.86 

7.41(  -8%) 

5.85(-15%) 

HgSe 

4.97 

5.77(+16%) 

5.11(  +3%) 

HgTe 

4.23 

4.64(  +6%) 

4.49(  -1%) 

Table  4.  Comparison  of  Cu— Ci2,  in  units  of  1011ergs/cm3,  for  various  zincblende 
semiconductors  using  the  local  density  (LD)  approximation  and  the  gradient  correc¬ 
tion  (GC)  to  the  local  density.  Shown  for  comparison  are  the  experimental  values. 
Values  shown  in  parentheses  are  the  percentage  difference  from  the  experimental  val¬ 
ues. 


Compound 

Experiment 

LD 

GC 

Si 

10.27 

10.26(  0%) 

9.83(  -4%) 

Ge 

8.19 

7.84(  -4%) 

7.47(  -9%) 

AlSb 

4.43 

4.64(  +5%) 

4.42(  0%) 

GaP 

7.87 

8.14(  +3%) 

7.70(  -2%) 

GaAs 

6.55 

6.42(  -2%) 

6.07(  -7%) 

GaSb 

4.95 

4.69(  -5%) 

4.46(— 10%) 

InP 

4.46 

4.63(  +4%) 

4.45(  +0%) 

InAs 

3.80 

3.86(  +2%) 

3.68(  -3%) 

InSb 

4.82 

2.99(—38%) 

2.96(— 39%) 

ZnS 

3.54 

4.40(+24%) 

4.40(+24%) 

ZnSe 

3.56 

3.78(  +6%) 

3.80(  +7%) 

ZnTe 

3.06 

3.67(+20%) 

3.42(-fl2%) 

CdTe 

1.68 

2.06(+23%) 

2.01(+20%) 

HgS 

1.91 

1.63(— 15%) 

1.65( — 14%) 

HgSe 

1.70 

1.49(— 12%) 

1.62(  -5%) 

HgTe 

1.82 

1.80(  -1%) 

1.81(  -1%) 
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ZnSe  — ►  VznSe  -f  Zn jTtt  atom) 

it  is  this  feature  of  the  GC  in  which  we  are  most  interested.  It  appears  that  the  GC 
is  most  important  not  in  the  calculation  of  the  energy  of  the  solid,  but  rather  in  the 
energy  of  the  free  atom.  Our  work  on  other  elemental  solids  in  the  periodic  table  -  for 
example  solid  antimony,  aluminum,  zinc,  and  tellurium  -  will  allow  us  to  substantiate 
this  observation.  The  work  on  the  gradient  correction  has  also  been  supported  by 
ONR  contract  N00014-91-C-0139.  Although  the  majority  of  the  computations  we  plan 
to  do  to  test  the  gradient  correction  have  been  completed,  we  expect  to  complete  our 
analysis  and  interpretation  of  the  results  in  the  next  quarter. 

1.2  DEGENERACY  FACTORS  FOR  DONOR  AND  ACCEPTOR  LEVELS 

One  of  the  ingredients  that  enters  into  the  calculation  of  the  defect  densities  is 
the  degeneracy  of  the  particular  ionization  state  of  the  defect.  While  the  formulas  for 
these  degeneracies  are  well  established  for  the  shallow-level  impurities  with  a  single 
ionization  state,  the  degeneracies  of  the  deep  state,  such  as  those  that  arise  from  the 
localized  potentials  of  the  native  defects,  have  not  been  worked  out.  We  have  com¬ 
pleted  an  analysis  of  the  degeneracies  for  the  states  for  the  various  donor  and  acceptor 
levels  of  the  vacancies,  antisites,  and  interstitials  in  the  zincblende  semiconductors, 
and  incorporated  our  results  into  our  codes  for  calculating  the  defect  densities.  The 
details  of  this  analysis  will  be  included  in  a  long  paper  under  preparation  for  publi¬ 
cation. 

1.3  TELLURIUM  INCLUSION  IN  MCT 

We  have  begun  to  investigate  the  formation  of  tellurium  inclusions  and  their  re¬ 
lationship  to  the  mercury  vacancy  and  tellurium  antisite  densities  in  HgCdTe.  In  the 
analysis,  we  have  used  our  results  on  the  defect  concentrations  in  HgCdTe  as  well  as 
some  of  the  literature  on  self-diffusion  in  HgCdTe.  Although  this  problem  necessitates 
more  information  about  diffusion  than  is  currently  available,  we  have  used  available 
information  to  construct  a  description  of  the  inclusions  and  their  environments.  Our 
tentative  conclusions  are 

•  Tellurium  antisites  diffuse  via  a  mercury  vacancy  mediated  process, 
and  therefore  the  diffusion  coefficient  is  proportional  to  the  product  of 
the  mercury  vacancy  and  tellurium  antisite  concentrations. 

•  The  migration  energy  of  the  tellurium  antisite  will  be  comparable  to 
that  of  the  mercury  vacancy,  while  the  overall  diffusion  coefficient  will 
be  much  smaller. 

•  In  HgCdTe,  tellurium  diffuses  via  tellurium  antisites  only;  tellurium 
interstitial  and  tellurium  vacancy  diffusion  is  negligible  because  of  the 


negligible  populations  of  these  species.  This  conclusion  is  strongly  sup¬ 
ported  by  our  work  on  the  native  defects  in  HgCdTe. 

•  A  complex  consisting  of  two  tellurium  antisites  (each  a  single  donor) 
and  one  mercury  vacancy  (a  double  acceptor)  may  be  the  nucleation 
site  of  a  tellurium  inclusion.  This  complex  may  be  well  bound  because 
of  its  charge  neutrality. 

•  There  may  be  an  enhanced  diffusion  of  tellurium  towards  a  tellurium 
inclusion  due  to  a  coordination  of  Coulomb  and  strain  terms. 

•  If  the  tellurium  antisite  diffusion  is  very  slow,  as  we  expect,  the  tel¬ 
lurium  antisites  may  not  equilibrate  during  the  low-temperature  mercury- 
saturated  anneals,  and  may  in  fact  be  the  residual  donor  that  becomes 
apparent  when  the  mercury  vacancy  concentration  is  reduced  below  the 
level  of  1018  cm-3. 

•  We  expect  that  the  tellurium  antisites  are  stripped  out  of  the  material 
in  the  vicinity  of  a  tellurium  inclusion  by  the  inclusion  formation  pro¬ 
cess.  Thus,  during  a  low-temperature  anneal,  the  material  adjacent 
to  the  inclusion  will  remain  mercury  vacancy  doped  p-type  after  the 
low-temperature  anneal,  while  the  bulk  will  turn  n-type  because  of  the 
relatively  high  concentrations  of  tellurium  antisites  there. 

•  If  these  p-type  regions  about  inclusions  intersect  one  another  to  create 
p-type  percolating  paths,  this  may  account  for  the  anomalously  low 
mobility  in  LPE  material. 

While  some  of  the  above  conclusions  are  tentative,  the  picture  they  present  is  con¬ 
sistent  with  our  theoretical  results  and  the  diffusion  experiments.  Much  more  work 
needs  to  be  done  to  complete  the  picture.  In  the  short  term,  examining  the  binding 
energy  of  the  tellurium  antisite  -  mercury  vacancy  complex  will  give  us  some  further 
insight  into  the  diffusion  of  the  tellurium  antisite. 

2.  DISLOCATIONS  IN  MCT 

We  completed  a  paper  on  our  work  on  dislocations  in  HgCdTe  and  the  way  in 
which  they  impact  device  performance.  A  copy  of  the  paper,  which  will  be  presented 
at  the  1993  Materials  IRIS  meeting,  is  included  as  Appendix  A. 

3.  TRANSPORT  PROPERTIES  OF  HgCdTe 

We  have  calculated  the  absorption  coefficient,  Fermi  energy,  Hall  coefficient,  and 
electron  mobility  with  Fermi- Dirac  statistics  (FD)  and  accurate  pseudo  potential 
band  structures  fine  tuned  with  tight-binding  (TB)  corrections.  We  have  studied  the 
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effect  of  the  Maxwell-Boltzmann  (MB)  and  parabolic  band  approximation  on  elec¬ 
tron  transport  coefficients  and  on  ways  in  which  physical  parameters  are  extracted 
from  experiments.  The  values  interpreted  from  experiments  depend  crucially  on  var¬ 
ious  approximations  such  as  effective  mass,  Maxwell-Boltzmann  statistics,  and  the 
collision  time.  Our  main  conclusions  to  date  are 

•  Approximating  the  Hall  factor  by  unity  over  a  wide  range  of  carrier 
concentrations  and  temperatures  is  accurate  for  most  applications.  An 
error  of  about  30%  is  expected  at  high  temperatures  and  low  carrier 
concentration  (1014  cm-3),  and  about  20%  is  expected  at  low  T  and 
high  carrier  concentration(H)1*  cm-3). 

•  The  variation  of  band  gap  with  temperature  is  nonlinear  and  is  faster 
at  low  temperature.  The  value  at  22%  mercury  concentration  is  about 
10  meV  smaller  than  those  usually  quoted. 

•  The  mobility  calculated  from  a  full  solution  to  the  Boltzmann  trans¬ 
port  equation  with  FD  statistics  can  explain  the  hump  near  40  K  and  is 
in  good  agreement  with  experiments.  Once  the  other  scattering  mech¬ 
anisms  such  as  acoustic  phonon  and  alloy  disorder  are  included,  the 
agreement  is  expected  to  be  much  better. 

This  work  will  be  reported  at  the  1993  Materials  IRIS  meeting.  A  copy  of  the  paper 
to  be  presented  at  that  meeting  is  included  as  Appendix  B. 

4.  WIDE-GAP  1 1- VI  COMPOUNDS  (ZnSe  AS  PROTOTYPE) 

Our  computational  resources  have  been  consumed  by  our  work  on  the  gradient 
correction  this  quarter.  As  discussed  above,  this  work  on  the  gradient  correction  was 
necessary  to  provide  the  foundation  for  further  work  on  defect  formation  energies. 
We  have  just  resumed  our  work  on  the  calculation  of  the  defect  formation  energies  in 
ZnSe.  In  the  last  quarterly,  we  reported  on  our  preliminary  calculation  of  the  defect 
formation  free  energies.  We  have  resumed  these  calculations  to  attempt  to  resolve 
discrepancies  with  other  published  work  and  to  continue  the  calculation  of  the  energy 
of  the  localized  states  associated  with  the  defects. 

5.  NONLINEAR  OPTICAL  MATERIALS  (LiNb03  AS  THE  PROTOTYPE) 

In  the  previous  quarter,  we  obtained  the  first-principles-based  band  structure 
for  LiNb03.  With  this  result  in  hand,  we  have  begun  our  investigation  of  the  de¬ 
fect  structure  of  this  compound.  The  difficulty  of  growing  stoichiometric  LiNb03 
is  well  known.  Most  growth  methods  yield  congruent  lithium  niobate  with  Li:Nb 
*  48.45:51.55.  Although  properties  such  as  lattice  parameters,  density,  dielectric 
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constant,  and  refractive  indices  depend  only  slightly  on  the  Li  to  Nb  ratio,  other  im¬ 
portant  properties  such  as  optical  absorption,  curie  point,  birefringence,  and  nuclear 
magnetic  resonance  (NMR)  line  width  are  altered  strongly  by  that  ratio. 

The  deviation  from  stoichiometry  will  certainly  introduce  defects.  Understanding 
the  defect  structure  of  this  compound  is  of  great  device  interest  and  the  focus  of  the 
present  study.  There  are  several,  often  contradicting,  models  of  the  defect  structure 
that  have  been  proposed  exclusively  from  experimental  information.  The  models 
include  Nb5+  residing  on  Li+  sites  with  four  lithium  vacancies  (for  charge  neutrality), 
arrangement  of  cations  along  the  polar  axis  out  of  normal  sequence  with  some  niobium 
atoms  on  lithium  sites,  excess  Nb  with  Nbu,  Li+  and  O-2  vacancies  with  no  stacking 
fault,  and  oxygen  vacancies  with  one  or  two  electrons  in  partially  reduced  lithium 
niobate.  All  these  models  explain  a  specific  experiment,  but  an  overall  defect  picture 
is  not  obtained.  Due  to  sensitive  dependence  of  various  defects  on  growth  conditions, 
various  characterization  techniques  and  possibly  small  activation  energy  difference,  a 
complete  defect  structure  is  still  not  available.  Once  such  a  knowledge  is  gained,  this 
material  can  be  engineered  more  efficiently  for  various  applications. 

We  have  scrutinized  the  models  proposed  in  the  literature  to  identify  the  defects 
in  LiNbOa  for  our  initial  study.  We  have  identified  the  lithium  vacancy  (Vu)>  the  nio¬ 
bium  antisite  (Nbu),  and  the  oxygen  vacancy  (Vo)  as  the  first  candidates  for  study, 
and  have  started  with  the  lithium  vacancy.  The  problem  is  approached  in  a  supercell 
approach  within  the  full-potential  linearized  muffin-tin  orbital  (LMTO)  method,  dis¬ 
cussed  extensively  in  the  previous  quarterlies  and  our  proposal.  We  have  started  with 
the  30-atom  unit  cell  containing  one  lithium  vacancy.  Although  this  corresponds  to 
6%  vacancy  concentration  we  have  found  in  our  work  on  other  semiconductors  that  the 
energies  from  these  calculations  are  nearly  converged  to  the  low  concentration  limit 
for  the  cell  size  quoted.  When  our  initial  survey  of  the  primary  defects  is  complete 
and  the  most  important  defects  identified  in  our  study,  we  will  refine  our  calculation 
for  the  most  important  defects.  The  structure  about  the  lithium  vacancy  is  relaxed 
to  minimize  the  defect  total  energy,  we  find  that  the  three  oxygen  atoms  in  the  plane 
of  removed  lithium  move  close  to  each  other.  We  have  encountered  some  numerical 
problems  in  calculating  this  relaxation  energy,  but  have  developed  some  strategies  to 
circumvent  the  problems. 


6.  WORK  PLANNED 

Four  papers  will  be  presented  in  the  next  quarter,  two  at  the  1993  Meeting  of  the 
IRIS  Specialty  Group  on  Infrared  Materials,  in  Bedford,  MA,  and  two  at  the  1993 
U.S.  Workshop  on  the  Physics  and  Chemistry  of  Mercury  Cadmium  Telluride  and 
other  m  materials  in  Seattle,  WA.  We  also  plan  to  submit  versions  of  these  papers  to 
the  Journal  of  Applied  Physic*.  Approval  for  publication  and  presentation  of  these 
papers  has  been  obtained.  We  will  be  finishing  our  work  on  the  gradient  corrections 
and  writing  a  paper  on  the  results;  this  paper  will  be  essential  to  the  support  of  our 
work  on  the  defects  in  HgCdTe,  ZnSe,  and  LiNbOa .  We  expect  to  complete  our  paper 
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on  defects  in  HgCdTe  for  submission  to  Physical  Review  B  in  the  next  quarter  as  well. 
We  will  begin  our  study  of  defects  in  CdTe  this  next  quarter  and  continue  our  work 
on  defects  in  ZnSe  (in  particular  the  localized  states)  and  LiNbOs- 
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ABSTRACT 

Dislocations  crossing  a  junction  in  HgCdTe  have  little  effect  on  detector 
responsivity,  but  are  known  to  reduce  the  zero  bias  impedance  RoA  and  increase 
the  leakage  current,  especially  at  low  temperatures  where  RoA  is  dominated  by 
tunneling  and  GR  processes1*2.  We  have  calculated  the  core  charge  and 
piezoelectric-induced  fields  associated  with  60*  dislocations  since  they  are 
superimposed  on  die  built-in  and  applied  junction  fields  that  control  the  currents. 
The  screening  of  the  dislocation  fields  in  the  neutral  region  is  nearly  complete  and 
is  consistent  with  die  dislocations’  small  effect  on  responsivity.  Their  impact  in 
the  space  charge  region  is  found  to  be  significant  and  consistent  with  the  nonlinear 
dependence  of  performance  on  dislocation  density,  as  characterized  by  RoA. 


1.  INTRODUCTION 


The  p-on-n  Hgi.xCdxTe  (MCT)  double-layer  heterojunction  (DLHJ)  focal  plane  arrays 
(FPAs)  shown  in  Figure  1  are  currently  leading  the  art  in  long-wave  infrared  (LWIR) 
applications.  The  most  effective  structures  to  date  are  those  grown  by  liquid-phase  epitaxy 
(LPE)  cm  lattice-matched  CdZnTe  substrates.  Typical  values  of  active  (n  layer)  and  cap  (p  layer) 
compositions  are  x  =  0.224  and  x  £  0.28,  respectively.  The  n  layers  are  typically  indium-doped 
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Figure  1 .  Cross-sectional  view  of  DLHJ  mesa  diode 

Source:  R.E.  DeWames  et  al.,  An  Assessment  of  HgCdT e  and  GaAs/GaAIAs  Technologies  for  LWIR 
Infrared  Imagers.  SPIE,  Voi.  1735,  p.  6, 1992. 

in  the  growth  process,  with  carrier  concentrations  as  low  as  permitted  by  series  resistance  and 
capacitance  concerns,  but  usually  realized  at  >  2  x  1015  cm*3.  The  cap  p+  layers  arc  typically 
arsenic-doped  with  p  =  2  x  1017  cm*3.  The  metallurgical  p-n  junction  is  positioned  on  the 
narrow-gap  side  of  the  heterojunction,  and  the  grading  width  is  chosen  to  minimize  the  barrier  in 
the  valence  band  that  occurs  in  such  structures3,  as  the  barrier  reduces  the  collection  efficiency  of 
photon-induced  minority  carriers  (holes)  in  the  n  layer  from  the  back-lit  array. 

There  is  ample  experimental  evidence  that  dislocations  degrade  the  performance  of  DLHJ 
detectors  and  may  be  the  defects  currently  limiting  die  capabilities  of  the  devices.  The 
measurements  that  relate  various  material  parameters  pertinent  to  detector  performance  to 
dislocation  density  have  been  reported  by  several  groups.1*2  The  important,  interrelated 
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parameters  are  minority  carrier  lifetime,  zero-bias  junction  impedance  (RoA  product),  and  1/f 
noise. 

Dislocations  are  thought  to  influence  the  electronic  properties  of  semiconductors 
primarily  through  their  electric  fields.  In  this  paper  we  report  the  results  of  analyses  of  the  space 
charge  and  piezoelectric  fields  associated  with  dislocations  in  an  attempt  to  understand  the 
mechanism  through  which  they  affect  the  responsivity  and  dark  currents  at  low  temperature. 


2.  SUMMARY  OF  EXPERIMENTAL  RESULTS 


In  the  woik  to  be  discussed,  dislocation  densities  are  measured  by  etch  pit  densities 
(EPDs).  The  etch  pits  on  the  film  surface  are  associated  with  threading  dislocations  that  are 
assumed  to  penetrate  the  films.  The  relation  of  these  dislocations  to  misfit  dislocations  and  loops 
formed  at  growth  interfaces  is  not  fully  understood  in  all  cases. 

Shin  et  al ,2  have  studied  the  relationship  between  dislocation  density  and  minority  carrier 
lifetime  in  epitaxial  MCT  films  grown  by  LPE,  metalorganic  chemical  vapor  deposition 
(MOCVD),  and  molecular  beam  epitaxy  (MBE).  At  small  dislocation  densities  (as  measured  by 
etch  pit  density)  the  minority  carrier  lifetime  is  not  affected  by  the  presence  of  dislocations,  and 
the  temperature  and  carrier  concentration  dependence  of  the  lifetime  suggests  Auger 
recombination  at  high  temperature  (above  150  K)  and  Shockley-Read  recombination  for 
T  <  150  K.  However,  even  in  samples  with  very  low  EPDs  (down  to  pixels  with  no  etch  pits), 
the  lifetime  and  the  RoA  product  exhibit  great  variability,  sometimes  exceeding  an  order  of 
magnitude.  In  these  high-quality  samples,  the  transition  from  Auger  to  trap-assisted 
recombination  can  occur  at  temperatures  as  low  as  60  K.4  We  take  this  observation  as  evidence 
that  recombination  centers  not  associated  with  threading  dislocations  may  be  the  defects 
imposing  the  ultimate  device  performance  limits.  In  the  low  temperature  range  of  interest  for  IR 
detectors,  for  EPD  densities  above  the  mid  10' 5  range,  minority  carrier  lifetime  decreases  with 
dislocation  density  at  an  EPD'2/3  to  EPD*3/4  rate.  Shin  et  al .2  note  that  the  observed  behavior  is 
consistent  with  dislocations  acting  as  Shockley-Read  recombination  centers  with  a  trap  level 
60  meV  above  the  valence  band  edge  and  a  trap  density  proportional  to  EPD'3/4.  Although  not 
explicitly  stated  in  this  paper,  the  results  reported  suggest  that  dislocations  active  as 
recombination  centers  are  likely  to  carry  a  core  charge.  Similar  results  have  been  reported  by 
Baianskii  et  al. 5  who  find  that  mobile  dislocations  introduced  in  bulk  MCT  by  indentation  give 
rise  to  a  trapping  level  at  Ev  +  50  meV. 
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Johnson  et  al .*  have  measured  the  dependence  of  RoA  in  arrays  of  photovoltaic  diodes  on 
dislocation  density  by  deliberately  using  plastic  deformation  to  introduce  dislocations  in 
localized  regions  and  carrying  out  the  measurement  on  both  high  and  low  EPD  diodes  of  the 
same  array.  They  find  that  the  resistance  is  relatively  insensitive  to  EPD  at  the  lowest  dislocation 
densities,  decreases  linearly  at  higher  densities,  and  decreases  quadratically  at  the  highest 
densities.  The  nonlinear  behavior  sets  in  at  small  dislocation  densities  (mid  10s  cm*2  or  less)  at 
low  temperature  (40  K)  and  at  higher  densities  (MO6  cm*2)  at  78  and  120  K.  Measurements  of 
the  temperature  dependence  of  the  leakage  current  of  diodes  with  different  dislocation  densities 
show  that  the  high  temperature  (diffusion  regime)  RoA  is  not  affected  by  the  presence  of 
dislocations,  but  that  as  the  temperature  is  reduced  below  ISO  K  the  RoA  values  diverge, 
differing  by  several  orders  of  magnitude  at  40  K  and  below.  The  authors  interpret  the  low 
temperature  leakage  current  as  generation-recombination  current  They  adopt  a  model  (Zolper 
and  Barnett6)  for  estimating  the  reduction  of  the  minority  carrier  diffusion  length  by 
recombination  at  dislocations,  and  then  account  for  the  inverse  quadratic  increase  of  the  short 
circuit  current  at  high  dislocation  densities  by  the  ad  hoc  assumption  that  closely  spaced 
dislocation  pairs  contribute  a  conductance  that  exceeds  the  sum  of  the  conductances  of  two 
individual  dislocations.  We  shall  point  out  below  that  such  nonlinear  behavior  may  be  a  feature 
of  tunneling  currents  associated  with  dislocations.  The  same  work  also  reports  on  measurements 
of  1/f  noise.  The  noise  current  was  found  to  be  proportional  to  the  leakage  current,  leading  to  the 
conclusion  that  the  dislocations  act  only  as  indirect  noise  sources  by  increasing  the  leakage 
current. 


3.  CALCULATION  OF  DISLOCATION  PROPERTIES 


We  wish  to  distinguish  between  dislocations  that  intersect  the  junction  and  those  that  do 
not  We  shall  refer  to  intersecting  dislocations  as  threading  dislocations.  Threading  dislocations 
eventually  intersect  the  surface  of  the  device,  and  therefore  are  observable  by  etch  pits.  Such 
dislocations  have  a  number  of  sources.  They  can  originate  in  the  substrate  and  propagate  through 
the  active  material  during  growth;  they  can  be  introduced  by  plastic  deformation  of  the  material; 
or  they  can  originate  from  lattice  or  thermal  mismatch  in  the  structure. 

3.1  DISLOCATION  CORES 

Classical  elasticity  theory  works  well  for  the  description  of  long-range  strain  fields  of 
dislocations  in  real  solids,  but  fails  when  it  conies  to  a  description  of  the  dislocation  core.  The 
description  of  dislocation  cores  necessarily  depends  on  the  microscopic  properties  of  the 
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materials.  If  we  view  the  60°  dislocation  as  the  removal  of  three  { 1 12}-type  half-planes  of 
atoms,  followed  by  the  displacement  and  rejoining  of  the  crystal  across  the  missing  plane,  the 
line  that  describes  the  termination  of  the  missing  plane  is  the  dislocation  line.  The  dislocation 
core  is  the  cylindrical  volume  centered  at  the  dislocation  line;  in  that  core,  bonds  are  broken, 
displacements  of  atoms  are  large,  and  displacements  cannot  be  described  by  linear  elastic  theory. 
Because  broken  bonds  are  present  in  the  cores,  localized  states  can  be  introduced  in  the  band  gap. 
Although  reconstruction  may  take  place  to  move  these  states  out  of  the  gap,  it  is  not  clear 
whether  this  occurs  in  MCT. 

For  the  undissociated  60°  dislocation  in  MCT,  the  local  states  of  dislocation  cores  will 
depend  on  the  plane  on  which  the  extra  half-planes  of  atoms  terminate  (i.e.,  whether  they  are 
cation  or  anion  terminated),  and  whether  the  dislocation  belongs  to  the  shuffle  or  glide  set. 
Dislocations  of  the  glide  set  in  the  group  IV  and  III-V  semiconductors  are  frequently  found  to  be 
dissociated  into  two  Shockley  partials  separated  by  a  stacking  fault.  One  study  on  GaAs  found 
that  as  many  as  80%  of  the  60°  dislocations  in  plastically  deformed  GaAs  were  dissociated7.  It  is 
not  clear  whether  the  majority  of  the  60°  dislocations  in  MCT  are  dissociated,  but  is  likely  that 
they  are.  Since  the  primary  effects  of  the  space  charge  and  piezoelectric  fields  considered  here 
occur  at  the  longer  range,  where  classical  elasticity  is  valid,  the  effects  considered  in  this  work 
differ  little  between  an  undissociated  dislocation  and  a  pair  of  partials. 

&2  CHARGED  DISLOCATIONS 

The  Coulomb  fields  associated  with  electrically  charged  dislocation  cores  can  act  as 
scattering  and  recombination  centers  and  may  attract  or  repel  charged  point  defects.  The 
presence  and  amount  of  charge  on  a  dislocation  core  depends  on  the  distribution  of  localized 
levels  in  the  bandgap  in  the  core  and  on  the  Fermi  level  in  the  semiconductor.  Once  the  core 
structure  of  the  dislocation  has  been  determined,  we  can  take  a  phenomenological  view  and 
consider  only  the  electric  field  and  potential  generated  by  a  given  core  charge.  The  effective 
radius  of  the  field  is  taken  to  be  that  distance  from  the  dislocation  core  within  which  the  field  is 
large  enough  to  move  a  minority  carrier  to  the  core  in  a  carrier  lifetime  as  measured  in  the 
dislocation-free  material.  In  this  way  we  define  a  volume  within  which  the  lifetime  and  mobility 
are  affected  by  the  presence  of  a  charged  dislocation.  Only  if  the  affected  volume  surrounding 
all  the  dislocations  is  a  significant  fraction  of  the  total  volume  can  we  expect  their  presence  to 
affect  lifetime  and  mobility  and  hence  RoA.  This  then  relates  the  depletion  radius  R  or  screening 
radius  t  to  a  threshold  of  the  dislocation  density  no- 

We  represent  the  dislocation  core  by  a  cylindrical  charge  density  -p  in  a  cylinder  of 
radius  a  in  n-type  material  with  donor  density  n.  The  potential  distribution  in  the  depletion 
approximation  is  given  by  the  Poisson  equation  with  charge  density  p: 
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p  =  (-  p  +  n)e 


r<a 


=  ne  a  <  r  <  R 

=  0  R  <  r 

where  R  is  the  depletion  radius,  related  to  a  by  the  charge  neutrality  requirement  pa2  =  nR2.  By 
integrating  the  Poisson  equation  we  find  for  the  potential 

<Kr)  =  0(a)  +  ^R2ln^  -  -j(r2  -  a2)j 

and  taking  the  potential  in  the  bulk  (r  £  R)  as  reference,  we  can  evaluate  the  potential  at  the  core 
radius 

*(a)  =  g[i(R2-a2)-R2lnf.]. 

Using  a  dislocation  core  radius  of  5  x  10*8  cm  and  assuming  the  bands  bend  by  the 
bandgap — that  is,  by  a  potential  (relative  to  the  bulk)  at  the  core  radius  of  0.1  V  for  8-  to  12-pm 
material,  we  have  calculated  R  as  a  function  of  carrier  density,  as  shown  in  Figure  2.  It  turns  out 
that  the  radius  is  close  enough  to  the  Debye  length  to  make  the  depletion  approximation 
questionable. 

The  exact  space  charge  equation  is 


liL  £l£ 

r  drrdr 


— r[sinh(<p  -tpp)  -  sinhtpp] 
L2 


where 


,2  EkT  _  Ef-Ej  ..^Nq-Na 
^  kT’ L  T“>  ^  — r^r~  smh 


c2n; 


kT 


n: 


For  9  small  this  is 


1  d  dtp 
r  drrdr 
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Figure  2.  Depletion  radius  R  and  screening  radius  l  (meters)  at  40, 78,  and  300  K  as 
functions  of  carrier  density  in  cnr3 


with  the  extrinsic  Debye  length 


■  2  L2  elcT 
"  e2N[) 


2.77x10 5  „  .  .3 

- r: - ,  Nq  in  cm  . 


Then  with  x  «  we  have  ^  ^  x  ^  =  <p,  the  zero-order  modified  Bessel  equation.  The  solution 

of  interest  is  the  decreasing  Hankel  function  Ko(x). 

The  depletion  approximation  should  be  valid  for  r  <  R  -  Lp.  This  does  not  leave  much 
room  since  R  and  Ld  are  comparable.  One  way  to  proceed  is  to  let  <p  =  CKo(x)  and  to  evaluate  C 
from  4>(a)  at  x  =  a/Lo-  Since  the  value  of  x  here  is  that  at  the  core  radius,  x  is  small  and  we  can 
use 


Ko(x)  =  ln(yx)  =  0.116-lnx, 


so 

Ko^).°.Il6  +  ln^).^ 
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to  determine  C.  The  effective  scattering  radius  of  the  dislocation  can  be  obtained  from  the  radius 
at  which  the  field  is  reduced  to  the  threshold  value  discussed  above,  here  denoted  E. 


E 


dr  e  dr  eLp  dx  v  £ 


where  n  is  the  carrier  density.  Here  the  argument  of  the  Bessel  function  corresponds  to  the 
screening  radius,  so  x  is  large,  and  we  can  use 


Finally  we  obtain  the  screening  radius  £  =  Ldx  from  the  solution  of 


The  threshold  field  E  is  now  determined,  as  discussed  above,  from  the  mobility  p  and  lifetime  x 
in  dislocation-free  material 


E  =  U\ix  *  Lox/px . 

This  is  substituted  in  the  immediately  preceding  expression,  and  the  resulting  expression  solved 
for  the  normalized  screening  radius  x.  The  screening  radius  f  as  a  function  of  n  is  plotted  for 
T  «  40  K,  78  K,  and  300  K  in  Figure  2.  In  the  calculation  we  have  used  the  exact  values  of  the 
Bessel  functions  rather  than  the  approximations  show  above,  and  the  following  values  of 
mobility  and  lifetime: 


T,K 

p,  cm2/Vs 

T,S 

40 

1.9  x  105 

1  x  10"6 

78 

9  x  104 

3x10*7 

300 

lxlO4 

lx  10*7 

As  expected,  £  turns  out  to  be  larger  than  the  depletion  radius,  confirming  the  need  for  the 
more  accurate  screening  theory.  The  figure  shows  the  results  of  both  calculations:  the  radius  R 
of  the  field  region  surrounding  the  positively  charged  core  as  a  function  of  carrier  density 
calculated  on  the  basis  of  the  depletion  approximation,  and  the  screening  radius  £  of  the  region 
inside  which  the  electric  field  due  to  the  charged  core  exceeds  the  critical  value.  The  threshold  is 
consistent  with  the  experimental  observations1*4  that  the  RoA  values  and  minority-carrier 
lifetimes  decrease  rapidly  with  the  dislocation  density  no  when  it  exceeds  106  cm*2.  The 


A-9 


calculation  of  the  screening  radius  £  is  more  reliable  than  the  depletion  approximation,  since  it 
takes  into  account  the  thermal  excitations  of  the  carriers.  With  the  parameters  we  use,  the 
calculation  of  the  screening  radius  over  five  orders  of  magnitude  of  the  carrier  density,  and  over 
the  temperature  range  from  40  to  300  K,  results  in  values  of  £  that  range  from  about  5  Ld  at 
n  =  10*0  cm*3  to  10  Ld  at  n  =  1015  cm*3,  resulting  in  an  overall  rate  of  decrease  of  £  with  carrier 
density  at  a  rate  slightly  slower  than  Ld'1^2. 

The  carrier  density  in  the  neutral  region  of  the  extrinsic  device  material  generally  equals 
or  exceeds  1014  cm*3.  The  screening  length  at  this  carrier  density,  even  at  room  temperature,  is 
less  than  5  pm,  and  it  is  only  about  3  pm  at  40  K.  Since  the  mobile  carrier  density  in  the  neutral 
regions  generally  exceeds  1014  cm*3,  we  can  conclude  that  in  the  presence  of  mobile  carriers, 
especially  at  cryogenic  temperatures,  the  core  charge  fields  have  a  limited  range.  This  is  also 
consistent  with  the  observations  of  Johnson  et  al .*  that  the  responsivity  of  p  on  n  heterojunction 
devices  (dependent  on  diffusion  in  the  neutral  region)  is  little  affected  even  by  high  dislocation 
densities. 

The  situation  is  less  clear  in  the  vicinity  of  the  metallurgical  junction  of  the  p-n  diode 
detector.  In  this  region,  the  carrier  density  is  reduced  to  the  intrinsic  value  nj  under  zero  bias,  and 
below  that  with  reverse  bias.  We  estimate  a  value  of  nj  about  5  x  101 1  cm*3  for  Hgo.8Cdo.2Te  at 
40  K,  with  a  decrease  of  about  a  factor  of  2  for  every  5  mV  of  reverse  bias  caused  by  the 
separation  of  the  quasi-Fermi  levels.  This  reduction  of  the  carrier  density  is  confined  to  a  small 
region,  and  an  assessment  of  its  quantitative  effect  will  require  detailed  computation. 
Nonetheless,  it  is  worth  observing  that  we  estimate  a  screening  radius  of  27  pm  for  this  carrier 
density.  This  is  consistent  with  the  experimental  finding4  that  RoA  at  this  temperature  is  already 
decreasing  with  dislocation  densities  in  the  10s  cm*2  range,  corresponding  to  this  order  of 
separation  between  dislocation  lines.  The  literature  values  of  nj  -  1.2  x  1014  cm*3  at  77  K  and 
ni  *  4.0  x  1016  cm*3  at  300  K7  correspond  to  progressively  shorter  screening  lengths,  which  is 
also  consistent  with  the  experimental  observation  that  the  decrease  in  RoA  with  dislocation 
density  has  a  higher  threshold  at  the  higher  temperature. 

The  transition  from  the  neutral  to  the  depleted  region  is  not  abrupt,  so  that  the  screening 
is  a  function  of  both  the  distance  from  the  dislocation  core  and  the  distance  from  the  junction 
plane.  For  a  dislocation  perpendicular  to  the  junction  plane,  the  geometry  has  cylindrical 
symmetry,  making  it  effectively  two-dimensional.  The  geometry  is  more  complicated  for  a 
dislocation  inclined  relative  to  the  junction  plane. 

Because  its  form  is  geometrically  simpler  and  more  readily  visualizable  in  this  context, 
we  have  discussed  screening  entirely  in  connection  with  the  effect  of  core  charges.  It  should  be 
kept  in  mind  that  similar,  although  quantitatively  different,  considerations  apply  to  the  screening 
of  the  piezoelec trically  induced  potentials,  to  be  discussed  next. 
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3.3  DISLOCATION  STRAIN  FIELDS  AND  PIEZOELECTRIC  POTENTIALS 


The  MCT  crystal  structure  lacks  inversion  symmetry  and  is  therefore  piezoelectric.  In  a 
piezoelectric  crystal  the  strain  field  around  a  dislocation  may  generate  an  electric  field. 
Conversely,  an  electric  field  applied  to  a  piezoelectric  crystal  can  produce  a  stress  that  may 
deform  the  crystal.  This  self-consistent  field  problem  has  been  treated  rigorously  for  a  straight 
dislocation  in  a  homogeneous  insulating  crystal  by  Saada8.  The  general  analysis  was  extended 
by  Faivre  and  Saada9  to  include  the  effects  of  screening  by  free  carriers.  Several  qualitatively 
significant  conclusions  can  be  drawn  from  this  work: 


1.  The  strains  are  only  slightly  affected  by  the  electric  field;  that  is,  the  inverse 
piezoelectric  effect  is  small,  and  the  electric  field  can  be  calculated  directly  without 
imposing  the  requirement  of  self-consistency. 

2.  The  electric  field  in  the  vicinity  of  the  dislocation  core  can  be  quite  large,  and  the 
dislocation  can  interact  strongly  with  charges  (point  defects  or  mobile  carriers)  in  its 
vicinity. 

3.  In  a  semiconductor,  the  dislocation  preferentially  traps  carriers  of  a  particular  sign. 

4.  The  estimated  screening  distance  is  of  the  order  of  several  Debye  lengths. 

Booyens  and  Vermaak10  have  applied  Saada’s  analysis  to  (1 10)/(1 1 1)  dislocations  in  zinc  blende 
ffl-V  semiconductors  and  Booyens  and  Basson11  to  CdTe  and  MCT.  They  do  not  explicitly 
include  the  effects  of  screening  as  analyzed  by  Faivre  and  Saada8  in  these  calculations.  As  a 
result,  they  find  that  an  uncharged  dislocation,  in  die  absence  of  screening,  generates  a  net  radial 
field.  They  then  suggest  that  the  dislocation  core  is  likely  to  acquire  a  charge  that  compensates 
for  the  radial  field  and  so  minimizes  the  total  energy  of  the  dislocation. 


We  have  extended  the  analysis  to  the  commonly  observed  60°  dislocations  in  intrinsic 
MCT.  The  dislocation  lies  along  [T01]  and  has  a  Burgers  vector  b  at  60°  to  the  dislocation  line 
4-  xi,  is  chosen  parallel  to  the  edge  component  of  the  Burgers  vector  b;  X2  is  chosen  as  the 
normal  to  the  slip  plane;  and  X3  is  chosen  parallel  to  £.  The  transformation  matrix  from  principal 
axis  coordinates  is 


f  1 
V2 
lV3 


1  ^ 
V2  , 

V3  J 


and  thus  the  transformation  of  the  piezoelectric  tensor  is 

d’ijk  *  TipTjqT|ndpqj. 
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In  cubic  symmetry  and  principal  axis  coordinates,  all  the  nonzero  components  of  the 
piezoelectric  tensor  are  equal:  di23=di32=d23l=d2l3=d312=d32l=sd-  The  nonzero  components  of 
the  transformed  tensor  are 

'  1  1  *  2  . 

d133  =  d313  =  d331  “  “dlll  ”  ~J7a 


d112  “  dl21  ~  d211  “  d233  ~  d323  “  d332  ~ 


'  _  , 
d222  ~  77?  ’ 


or,  using  two-index  (Voigt)  notation, 


(22  1  ^ 

— ^rd  0  -prd  0  0  -- =d 

V6  V6  V3 

(d')=  — Lrd  -~d  — =d  0  0  0. 

<3  <3  VJ 

0  0  0  — ==d  -2=d  0 

V  VJ  V6 

The  piezoelectrically  induced  bound  charge  density  is12 

1  3 

p(r)  =  -  d'ijkejk^  with  ejiy  e* 

£ 

where  e  is  the  isotropic  dielectric  constant,  ejk  is  the  jk  component  of  the  strain  tensor,  and 
summation  over  repeated  indexes  is  implied. 

For  a  screw  dislocation  in  linear  anisotropic  theory,  the  displacements  are 


u3=— -*-tan 
2n 


■'*(?} 


ui  *  U2  =  0, 

where  bfis  the  magnitude  of  the  screw  component  of  the  Burgers  vector  and 


rc’44c’55  1  1 

—J - ;  C44  =  C44  —  jH;  c’55  =  C44 -  ^H;  H  =  2044  +  c12  - c„. 
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Here  cy  and  c’ij  are  the  components  of  the  elastic  tensor  in  the  principal  axis  and  transformed 
coordinate  systems  respectively.  The  nonzero  components  of  the  strain  tensor  are  [using  ey  = 
(uy  +  ujj)] 

ei3  =  e3i  and  e23  =  e32. 

Using  Eqs.  1  and  2  shows  immediately  that  for  this  case  r  =  0.  Even  in  anisotropic  linear 
elasticity  (assuming  the  dielectric  constant  isotropic)  the  screw  component  of  the  60*  dislocation 
induces  no  piezoelectric  effect  Therefore,  pure  screw  (1 10)  dislocations  will  also  show  no 
piezoelectric  activity  even  if  they  are  dissociated.  We  now  proceed  to  consider  the  edge 
component,  which  we  will  simply  denote  b.  In  order  to  arrive  at  a  transparent  expression,  we 
will  make  the  isotropic  elastic  approximation.  This  will  not  affect  the  order  of  magnitude 
calculations  that  follow. 

For  an  edge  dislocation,  with  v  Poisson’s  ratio 


b  -i(  *2 

u'=af[Bn  Ur 


l  + _ _ 1. 

>  2(l-v)(x2+x|)  J* 


1— 2v 


u2  =  ^-|  -^-ln(xX)  + - 7  , 

2*[4(l-v)  4U-vXx2,+X2> 

u3  =  0. 


2  2 
Xi-X2 


The  nonzero  strain  components  are 


b  (1-2v)x|+(3-2v)xiX2 
2n  2(v-1)(xi+X2>2 


c22~ 

e12  = 


b  (1— 2v)x2~(1+2v)xjX2 
2n  2(v-l)(x^+x|)2 

b  xix2~ xi 

2n  2(v-1)(x21+x22)2  ' 


Frau  Eqs.  2  and  1, 


P^r)  “  ”  (dl  1  fil  1.1+  1^124+ d21 1C1  l,2+d222*22,2  ) 


A- 13 


*o|— 


so  we  need  the  derivatives 


b  (2v+ 1  )x  t  x  2-K2v-3)x:j  x2 

e,u"5T  (v-ikx2+i1)3 

b  (2v-l  )x2~6x  jx|-K3-2v)xi 
1U  2n  2(v-1)(xi+X2)3 

b  (2v-1)x2+6xix|-(2v+1)xi 
^  2n  2(v-1)(x?+X2)3 

_  b  X2-6X|X2+X| 

ClU  2n  2(v-l)(x2+x2)3 


giving 

p(xltx2) =  cn.i~2ei2,i~€nj.+2enj) 

bd  (2v-3)x^-2V2  (2v+  1  )x  jX2+30x2x2  +2^2  (3— 2v  )x3  xj~ (2v+7  )x{ 
2*V3e  2(v-1)(x2+x|)3 

When  this  expression  is  converted  to  polar  coordinates  in  the  xix2  plane,  it  becomes  somewhat 
simpler 

bd  /  , 

P  *  iitVle  (p>  +  p2+p3  +  p4  +  p5) 

with 

-  ^2v~3)  y4  _  (2v-3)  cos40-4cos20+  3 
Pl=2(v-l)r«“2(v-l)  8r2 

-2V2(2v-l)  xy3  -2V2(2v-l)  2sin20  -sin40 
Pz_  2(v-l)  r6  *  2(v-l)  gr2 

30  x2y2  30  l-cos40 
P3=  2(v-l)  r6  =  2(v-l)  8r2 
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2V2(3-2v)  x3y  2V2(3-2v)  sin40+2sin20 
P4  2(v-l)  r6  "  2(v-l)  gr2 

— (2v+7)  x4  -(2vf7)  cos40+4cos20+3 

Ps  “  2(v-l)  r6  “  2(v-l)  8r? 


- - ~[5cos40-V2  sin40+2(  l+v)cos20-V2  ( 1-2  v)sin20] 

4;tV3e(l-v)  r2 


—  [a„cos(n0)  4  o„sin(n0)] 

1  A 


which  defines  a  useful  constant  C,  and  four  parameters  depending  only  on  Poisson’s  ratio.  This 
simple  expression  reveals  the  charge  density  to  have  an  inverse  square  decay  and  quadrupole  and 
hexadecapole  terms  only.  Note  in  particular  that  there  is  no  monopole  term  and  hence  no  net 
charge  induced  on  the  dislocation  line.  This  rather  transparent  expression  has  not  been  revealed 
in  previous  work10*1 

To  obtain  the  piezoelectric  potential  we  need  to  solve  the  Poisson  equation 

with  this  charge  distribution  as  the  source.  In  two  dimensions,  the  formal  solution  of  the  Poisson 
equation  is 


O(r,0)  =  P  r’dr’ [  d0'  In—  f— 1  cosm(0-0')  p(r') 

*j„  Jo  [  ■■>  j 

=  -Y—  f  d0’[sin(m0)sin(m0')+cos(m0)cos(m0')]x 

k  *~4TnjQ 
m=l 

x  T  [a„cos(n8'}+b„sin(n8')]  f" T-f ~~  ] 
mu  •'0  <  \<>) 

where  r>  and  r<  are  the  greater  and  lesser,  respectively,  of  r  and  r\ 

In  analogy  with  the  solution  for  the  elastic  field,  the  integrand  diverges  at  the  origin  and 
is  also  not  normalizable.  To  deal  with  this  difficulty,  we  adopt  the  same  procedure  that  is  used  to 
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avoid  the  divergence  of  the  elastic  energy  of  a  dislocation  at  the  origin  in  continuum  theory.  We 
set  cutoffs  at  an  inner  radius,  Ri,  and  an  outer  radius,  R2-  The  cutoff  at  Ri  removes  the  core 
singularity  and  the  cutoff  at  R2  allows  p  to  be  expanded  in  multipoles.  Later,  we  will  identify 
Rjas  the  core  radius  and  we  will  let  R2  — >  We  now  have 

<Kr,0)  =  ^  ^  anCosCnQHbnSinCnS)^  J  0')n”1dr'  +  rnj^  (r')-"-^' 

=  C{  J  ta4cos(40>+b4sin(4e)]^  j  -  J 

+j  la2cos(29>+-b2sin(20)]^  1  -  ]  j 


forr<R2.  Forr>R2 


O(r,0)  =  c|  J  [a4Cos(46)+b4sin(46)]  -  (— ■)  ] 

+j  [a2cos(20)+b2Sin(20)]-“^~  j  -  (-y-)  jj 


which  is  a  multipole  potential.  The  two  solutions  are  matched  at  R2.  We  can  now  let  R2  ->  00 
and  neglect  the  terms  in  Ri,  which  amount  to  a  small  correction  near  the  core,  to  get  a  potential 
valid  for  r»Rj 

CKr.O)  =  c|  -g-  [a4cos(40)f b4sin(40)]  +  y  [a2cos(20>+b2sin(20)]| 


which  is  independent  of  r.  This  rather  surprising  behavior  is  a  consequence  of  the  inverse  square 
decay  of  the  charge  density.  In  this  connection,  we  can  pursue  the  analogy  with  the  well-known 
elastic  solution  for  an  edge  dislocation  by  pointing  out  that  while  the  elastic  energy  per  unit 
length  is  given  by  (11  is  the  isotropic  shear  modulus) 


Hb2 

4n(l-v) 


In  ^2.-0. 
Ri 


\ 

5 

j 


the  piezoelectric  energy  per  unit  length  is  (where  C  was  defined  earlier,  and  taking  v  =  0.3) 


^C2 

L6 
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which  is  six  orders  of  magnitude  smaller — justifying  our  neglect  of  the  reverse  piezoelectric 
effect  In  both  cases,  the  first  term  comes  from  the  energy  in  an  infinite  cylinder  integrated  to  R2 
and  which  diverges  logarithmically,  and  the  second  term  is  the  effect  of  the  cylindrical  surface 
bounding  this  dislocation  in  a  finite  crystal  (the  image  term). 

For  a  good  numerical  calculation,  we  would  need  the  piezoelectric  constant  of  the  alloy, 
which  has  not  been  determined.  For  an  order-of-magnitude  estimate,  we  can  use  the 
piezoelectric  constant  of  CdTe,  3  x  10"6  C  cm"2.  The  Burgers  vector  is  4  x  10*8cm,  and 
Poisson’s  ratio  can  be  taken  as  0.3.  Using  these  values  we  have  computed  the  charge,  potential, 
and  field,  shown  in  Figure  3. 

It  should  be  understood  that  because  of  the  symmetry  of  the  strain  field  associated  with  a 
straight  dislocation,  several  of  the  possible  sources  of  the  piezoelectric  potential  are  identically 
zero.  In  the  reduced  symmetry  of  curved  dislocations  (e.g.,  misfits),  kinks,  and  jogs,  these 
components  of  the  potential  do  not  vanish  and  may  be  significant 


3.4  AVERAGING  NONLINEAR  EFFECTS 

The  results  we  have  obtained  afford  a  qualitative  insight  into  the  effect  of  the 
piezoelectric  potential  on  the  characteristics  of  a  pn  junction  threaded  by  a  dislocation.  The 
angular  variation  of  the  potential  will  alternately  add  to  and  subtract  from  the  built-in  potential  of 
the  junction,  so  that  the  barrier  seen  by  a  carrier  crossing  the  junction  will  depend  on  its 
trajectory.  Because  of  the  strongly  nonlinear  dependence  of  both  diffusion  and  tunneling 
currents  on  the  barrier,  the  effect  of  such  a  barrier  variation  does  not  average  to  zero  and  the 
currents  can  be  enhanced  by  a  large  factor.  To  illustrate  the  principle,  we  choose  as  an  example 
the  average  tunneling  probability  for  carriers  crossing  a  narrow  semiconductor  junction. 

The  band'to-band  tunneling  probability  has  the  form 

p  =  exp(-CEg3/2) 

where  Eg  is  the  bandgap  and  C  is  a  constant  containing  all  other  material  parameters. 

We  crudely  represent  the  effect  of  the  dislocation  potential,  which  makes  the  tunneling 
barrier  a  function  of  the  carrier’s  trajectory,  as  an  angular  variation  of  the  bandgap  of  the  form 
EgsEgo(  1  +Asin0)  (it  will  be  apparent  presently  that  the  detailed  form  of  the  angular  variation  is 
irrelevant  for  this  argument).  It  is  likely  that  this  representation  of  the  piezoelectric  potential 
introduced  by  a  dislocation  can  probably  be  justified  by  a  WKB  calculation. 
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(a)  Piezoelectric  charge  normal  to  the  dislocation  line.  Solid  contours  are  1 , 0.5 

0.1,  and  0.05 10 12  e/cm3  going  out  radially;  dotted  contours  are  the  negative 
of  these. 

tin] 


(b)  Ouadrupole  component  of  the  piezoelectric  potential. 


Figure  3.  Piezoelectric  properties  of  a  60°  dislocation  in  Hg0  8Cdo  2Te 
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(d)  Angular  variation  of  the  piezoelectric  charge  and  potential  normal  to 
the  dislocation  fine 

Figure  3  (continued).  Piezoelectric  properties  of  a  60°  dislocation  in  Hg0  8Cdo  2Te 
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If  A  is  small  compared  with  unity  we  have 


E^2-E^(l+|Asine) 

(this  approximation  is  taken  purely  to  simplify  the  mathematics  and  has  no  effect  on  the 
argument).  We  can  now  compute  an  average  of  p 

<P>=4Jq  exp[  -CE^(  1+|  AsinQ  )]  d0=|£  Jo  exp[-|cE^ Asinsj  d0=poIO  (§CE^?a) 

where  Iq(x)  is  the  modified  zero-order  Bessel  function  and  po  is  the  tunneling  probability  in  the 
absence  of  any  dislocations.  From  this  last  form  it  can  be  seen  that  the  average  tunneling 
probability  could  have  been  written  as 

<P>  =  PoIo(,|AInIPol)- 

It  can  be  seen  that  although  the  fractional  bandgap  “angular  modulation”  A  is  small,  the  argument 
of  the  Bessel  function  need  not  be,  and  the  factor  multiplying  the  unmodulated  bandgap 
transition  probability  can  be  quite  large.  It  is  also  evident  that  any  other  angular  variation  would 
lead  to  very  similar  results.  The  angular  variation  A  can  be  estimated  from  the  angular  potential 
plots  of  Figure  3  to  be  roughly  0.05  to  0.1.  The  parameters  determining  po  (e.g.,  effective  mass, 
doping,  temperature)  cover  a  wide  range.  For  typical  values  of  these  parameters  the  multiplying 
Bessel  function  ranges  from  near  unity  to  about  10,  so  <p>/Po  is  in  rough  agreement  with  the 
observed1  effect  of  a  single  dislocation  on  RoA  of  a  pixel. 

As  another  consequence  of  the  type  of  nonlinearity  examined  here,  the  effect  of 
overlapping  piezoelectric  potentials  of  neighboring  dislocations  is  expected  to  be  enhanced  over 
what  might  be  expected  from  a  linear  superposition.  This  may  be  a  possible  clue  to  the  observed 
inverse  faster  than  linear  falloff  of  RqA  with  dislocation  density. 


4.  CONCLUSIONS 


We  have  evaluated  the  electric  fields  associated  with  dislocation  charges  and  strain  fields, 
and  estimated  their  effects  on  junction  performance.  We  conclude  that  the  fields  associated  with 
core  charges  alone  are  not  likely  to  account  for  the  observed  degradation  of  junction 
performance,  particularly  when  the  lower  carrier  concentration  is  on  the  n  side  of  the  junctions  so 
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the  principal  minority  carriers  are  holes.  The  longer  range  piezoelectric  potentials,  which  we 
have  calculated  explicidy  for  the  common  60°  dislocation,  are  likely  to  have  a  larger  effect, 
especially  on  RoA  through  the  tunneling  leakage  current.  The  nonlinear  dependence  of  the 
junction  characteristics  on  dislocation  density  is  examined,  and  is  proposed  as  a  natural 
explanation  for  the  observed  nonlinear  behavior  of  RoA. 
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ABSTRACT 

Transport-related  properties  such  as  electron  mobility.  Hall  coefficient,  Fermi 
level,  and  energy  gap  are  calculated  with  accurate  analytical  band  structures  and 
Fermi-Dirac  statistics.  Calculated  values  differ  substantially  from  the  ones 
obtained  with  effective  mass  or  parabolic  band  structure  approximation  of 
Hgo.22Cdo.78Te  alloy. 


1.  INTRODUCTION 


The  experimental  results  on  electron  transport  properties  of  semiconductors  are  often 
compared  to  theory  that  rests  on  three  approximations,  namely,  parabolic  band  structures  for 
those  states  occupied  in  the  measurement,  Maxwell-Boltzmann  (MB)  statistics,  and  collision 
time  approximation  to  a  full  Boltzmann  gain-loss  equation.  These  approximations  are  made  to 
all  scattering  mechanisms,  whether  they  are  elastic  or  inelastic.  It  is  well  known  that,  even  in 
large-gap  materials,  die  constant  effective  mass  approximation  is  valid  only  very  near  (*  Eg/10) 
to  the  band  edge.1*2  This  approximation  has  been  recognized  to  be  particularly  poor  for  narrow- 
gap  materials  and  nonparabolic  corrections  calculated  in  the  k*p  formalism.3-5  This  correction 
is  substantial,  but  still  differs  considerably  from  our  more  accurately  calculated  band  structures. 
In  addition,  our  fit  of  the  conduction  band  to  an  analytical  function  makes  many  results 
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Note:  This  presentation  is  unclassified  in  its  entirely. 


transparent  and  simplifies  the  calculations.  As  the  Fermi  energy  can  easily  move  into  the 
conduction  band  of  lightly  doped  small-gap  materials,  the  form  of  the  Boltzmann  equation  with 
Fermi-Dirac  (instead  of  the  usual  MB)  statistics  must  be  used  to  obtain  accurate  transport 
coefficients. 

We  report  results  from  our  study  of  absorption  coefficient,  Fermi  energy.  Hall  coefficient, 
and  electron  mobility  calculated  with  Fermi-Dirac  statistics  (FD),  and  an  accurate 
pseudopotential  band  structure  fine-tuned  with  tight-binding  (TB)  corrections. 


2.  BAND  STRUCTURE 


Quantitatively  accurate  band  structures  of  group  IV  elements,6  III-V  compounds,7  and 
II- VI  compounds7*8  can  be  obtained  using  a  minimum  set  of  sp3  orbitals  in  semiempirical 
calculations.  First,  empirical  pseudopotential  form  factors  are  used  to  calculate  a  TB 
Hamiltonian,  H  in  the  minimum  set.  H  is  then  transformed  into  a  zeroth  order  Ho  in  an 
orthonormal  basis.  Then,  a  perturbative  Hamiltonian  having  a  first-neighbor  TB  form  is  added  to 
Ho  to  fine-tune  the  band  structure.  Because  long-range  interactions  are  included  in  this 
Hamiltonian,  the  measured  band  curvatures  are  correctly  reproduced.  This  procedure  is  followed 
for  both  HgTe  and  CdTe,  and  then  the  alloy  band  structures  are  calculated  in  the  coherent 
potential  approximation. 

We  focus  Hgo.22Cdo.78Tc  alloy  with  100-meV  band  gap  for  the  studies  reported  here. 

We  find  that  the  calculated  conduction  band  is  replicated  very  well  by  a  hyperbola, 

Ek  =  (yk2  +  c2)1/2-c  (1) 

where  y  and  c  are  adjusted  to  fit  the  calculated  band  structure  in  the  energy  range  of  interest. 
When  y  and  c  are  treated  as  constants  related  to  the  band  gap,  Eg  and  the  effective  mass,  this 
expression  reduces  to  the  same  nonparabolic  correction  form  obtained  in  the  k*p  method.4 
However,  the  numerical  value  of  y  and  c  are  not  same  as  ours.  For  example,  in  the  chosen  case,  y 
and  c  are  48.3  and  0.058  respectively,  whereas  the  corresponding  k*p  values  are  36.0  and  0.05. 
The  differences  are  found  to  be  large  enough  to  cause  a  noticeable  change  ir  the  band  structure. 
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The  band  structure  calculated  by  diagonalizing  the  Hamiltonian  is  shown  in  Figure  1  (thick,  solid 
lines).  We  can  see  that  the  fitted  band  structure  (thin  line)  agrees  quite  well  up  to  an  energy  of 
0.5  eV  from  the  conduction  band  edge.  Without  loss  of  accuracy,  in  the  studies  considered  here, 
Eq.  1  is  used  as  the  energy-dispersion  relation  in  the  transport  expressions.  Also  shown  in 
Figure  1  is  the  poor  reproduction  of  the  conduction  band  obtained  with  an  effective  mass 
approximation  (dashed  line). 

Two  qualitative  features  of  the  band  structure  in  Figure  1  that  impact  transport  properties 
should  be  noted.  First,  for  energies  E-Ec  greater  than  50  meV  where  the  shape  of  the  conduction 
band  is  nearly  linear  in  k,  the  group  velocity  is  a  constant  independent  of  the  k.  Then,  the 
density  of  states  (DOS)  increases  proportional  to  E  rather  than  E1/2  as  in  the  case  of  parabolic 
bands.  Clearly,  these  features  modify  the  transport  properties  of  electrons  occupying  these  states. 
As  we  will  show  in  the  following  section,  at  the  carrier  concentration  and  temperatures  often 
found  in  device  structures,  the  Fermi  level  falls  into  the  region  where  these  features  contribute. 


3.  FERMI  LEVEL 


The  calculation  of  Fermi  level  Ef  as  a  function  temperature  T  and  doping  concentration 
no  is  required  for  all  transport  calculations.  A  knowledge  of  temperature-dependent  gap  EgfT)  is 
essential  to  obtain  £p.  Ideally,  temperature  dependence  should  be  included  in  the  Hamiltonian 
from  which  the  variation  of  Eg  with  T  could  be  obtained.  We  have  developed  a  general  method 
to  incorporate  phonon  and  alloy  effects  into  the  same  CPA  formalism,9  but  such  an  approach  is 
not  attempted  here.  Instead,  we  use  the  expression10  given  by 

Eg  =  0.0954  +  0.0327  T/1000  (2) 

After  studying  the  effects  of  various  approximations  and  obtaining  trends,  the  calculations  is 
repeated  with  the  temperature-dependent  Hamiltonian. 

Thai,  the  ep  is  calculated  from  the  condition1 1  that  at  a  given  T  the  number  of  electrons 
in  the  conduction  band  is  the  sum  of  electrons  excited  from  the  valence  band  and  donor  levels. 
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In  this  study  where  the  modifications  caused  by  the  band  structures  are  being  emphasized,  the 
donor  states  are  assumed  to  be  located  at  the  bottom  of  the  conduction  band.  The  valence  and 
conduction  band  DOS  are  calculated  from  our  band  structure.  The  valence  band  DOS  yields  a 
hole  effective  mass  of  0.65.  The  £f  (measured  from  valence  band  edge)  as  a  function  of  T  and 
nj>  calculated  from  parabolic  and  hyperbolic  band  structures  are  given  in  Figure  2,  a  and  b 
respectively.  As  expected,  a  hyperbolic  band  makes  a  substantial  change  from  parabolic-band¬ 
generated  values.  The  impact  on  the  transport  properties  is  large  when  ep  is  located  at  energy 
where  the  band  changes  its  character  from  parabolic  to  linear.  In  addition,  when  ef  is  near  to  or 
greater  than  Eg,  the  absorption  cutoff  wavelength  is  strongly  influenced.  A  proper  account  of  the 
Moss-Bemstein  shift  must  include  these  accurate  band  structures.  We  also  note  that  the  Fermi 
level  and  the  measured  band  gap  affect  each  other.  Hence,  a  proper  data  reduction  of  cutoff 
wavelength  determination  of  the  energy  gap  has  to  be  done  self-consistently. 


4.  ENERGY  GAP 


To  emphasize  the  point  that  the  measured  value  of  Eg  and  its  T  dependence  is  sensitively 
dependent  on  the  knowledge  of  Ef,  we  carried  out  a  calculation  of  absorption  coefficient  a.  The 
absorption  coefficient  a  is  proportional  to  k2  pc  (1  -  fc)  pv  fv.  where  f  is  the  FD  distribution 
function,  p  is  the  DOS,  and  subscripts  c,  v  represent  conduction  and  valence  band,  respectively. 
The  k2  factor  arises  from  the  matrix  elements  that  are  in  the  expression  for  ef.  The  temperature- 
and  energy-independent  proportionality  constant  is  adjusted  to  agree  with  the  experimental  curve 
at  80  K  in  the  vicinity  of  a  given  value  of  a.  Nearly  linear  dependence  of  ln(a)  on  photon  energy 
was  observed12  for  T  between  SO  and  300  K  and  Hg  concentrations  near  0.22. 

Because  the  band  edges  are  broadened  by  impurity  and  phonon  scattering  (Urbach  tail  or 
broadening),  the  gap  cannot  simply  be  assigned  to  the  energy  wave  length  where  these  curves 
project  to  zero.  Hie  procedure  used  by  many  authors  is  to  assign  Xco  to  be  the  place  where  a  is 
500  cm-1  or  1000  cm-1.  In  the  spirit  of  this  procedure,  we  adjust  our  proportionality  constant  to 
fit  the  bi(a)  versus  E  curve  at  80  K.  As  seen  from  Figure  3a,  one  constant  fits  the  entire  curve. 
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Figure  2a.  Fermi  energy  as  a  function  ol  temperature:  parabolic  bands 
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Figure  3b.  Energy  gap  as  a  function  of  temperature 
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However,  at  higher  temperatures  the  tail  sets.  We  then  recalculate  the  band  structures  and  a  for 
different  band  gaps  (with  coiresponding  curvature)  until  the  curves  at  a  =  500  cm-1  for  all 
temperatures,  as  shown  in  Figure  3a.  The  energy  gap,  variation  with  temperature,  deduced  by 
following  this  procedure  is  shown  in  Figure  3b  (a  dashed  line).  Also  shown  in  Figure  3b  (solid 
lines)  for  comparison  are  effective  Eg  versus  T  variations  if  values  are  taken  from  the  data 
corresponding  to  a  500  cm-1  and  1000  cm-1.  Finally,  because  our  procedure  is  ad  hoc,  we 
repeated  our  fits  by  adjusting  to  a  =  100  cm-1  values;  the  calculated  energy  gaps  are  also  plotted 
in  Figure  3b  (a  dashed  line).  We  expect  that  the  proper  Eg(T)  curve  lies  somewhere  between  our 
two  fit  curves.  Note  that  both  curves  are  below  the  traditional  curves  and  the  variation  in  T  is  no 
longer  linear,  emphasizing  the  need  for  proper  calculation  of  Eg(T),  a(T)  so  that  more  reliable 
device  parameters  can  be  established. 


5.  HALL  COEFFICIENT 


The  carrier  density,  n  in  n-type  material,  is  normally  deduced  from  measurements  of  the 
Hall  coefficient  Rh,  given  by  re/en,  by  assuming  the  Hall  factor  re  is  unity.  If  one  uses  parabolic 
approximation  and  MB  statistics,  re  is  approximately  unity.  We  set  out  to  examine  the  effect  of 
removing  these  approximations,  using  the  correct  band  structures.  This  requires  generalizing  the 
Boltzmann  transport  equation  (BTE)  to  include  the  FD  distribution  function.  We  have 

=  £  [w(k,  k’)f(k')(l-f(k))  -  w(k\  k)f(k)(l  -  f(k'))]  (3) 

The  first  term  of  the  right  side  of  Eq.  3  is  the  gain  term,  and  the  second  one  is  the  loss  term.  In 
equilibrium,  the  left  side  of  Eq.  3  is  identically  zero,  and  f  becomes  the  equilibrium  FD 
distribution  function  fo  given  by 

fo(Ek)  =  (cP(Ek"CF)+l)~1 

where  f}  is  (kfiT)-1.  In  the  presence  of  the  electric  and  magnetic  field 


(4) 


ffll=w)  +  vf(k).l(E  +  VXB) 


(5) 


In  steady  state,  the  df(k)/dt  in  Eq.  5  vanishes.  In  the  small-field  regime,  we  can  linearize  f  and 
write  it  as  a  sum  of  fo  and  a  perturbation  fi(k).  Disregarding  the  derivative  of  fi(k)  and  after 
some  algebraic  manipulation,  Eq.  5  reduces  to 

Vf0(k)  (E  +  v  x  B)  =  Y  [W(k,  k')fi(k’)  -  W(k\  k)f(k)fi(k)]  (6) 

n  p 

where  the  renormalized  W  and  the  usual  transition  probability  per  unit  time  w  are  related  by 

W(k,k')  =  w(k,k’)[11~W^  (7) 

Note  that  for  elastic  scattering  W  and  w  are  equal.  However,  for  inelastic  cases  the  effect 
depends  on  whether  energies  at  k  and  k'  are  larger  or  smaller  than  ep.  If  both  initial  and  final 
energies  are  larger  (or  smaller)  than  ef,  only  small  correction  to  w  is  expected.  However,  if  the 
initial  state  is  above  ef  and  the  final  state  is  below  ef,  that  scattering  is  suppressed. 

In  the  collision  time  approximation,  the  gain  term  in  Eq.  7  is  neglected  and  the  effective 
collision  time  is 

(<)"*  =  £  w(k\k)  (8) 


Using  this  collision  time  approximation,  it  is  straightforward  to  obtain  the  expression  for  re: 


£k  fo  Ik  k2  ^  fckF  W  -  fo) 

.  (Ikk^tx^fod-fo))2  . 


(9) 


where  VfcEfe  is  ft  k.  The  calculated  re  with  parabolic  and  hyperbolic  band  structures  is  shown 
for  various  T  and  doping  densities  in  Figure  4,  a  and  b  respectively.  In  both  cases,  the  re  is 
approximately  1  for  higher  temperatures.  At  lower  temperatures,  it  shows  considerable  structure 
with  a  maximum  value  of  about  2  obtained  in  parabolic  approximation.  However,  when  more 
accurate  hyperbolic  band  structure  is  used,  the  variation  with  T  is  reduced  to  a  maximum  value 
of  about  12,  suggesting  that  reported  intrinsic  densities  may  be  smaller  by  approximately  20%. 
The  effect  of  collision  time  approximation  and  T  dependence  of  band  gap  on  the  values  of  re  still 
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Figure  4a.  Halt  factor  as  a  function  of  temperature:  hyperbolic  bands 
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need  to  be  studied  to  extract  correct  carrier  densities  from  Hall  measurements.  However,  our 
tentative  conclusion  is  that  approximating  re  to  be  1  is  better  than  we  had  expected. 


6.  DRIFT  MOBILITY 


The  formalism  developed  above  to  find  a  solution  to  BTE  with  FD  statistics  can  be  used 
to  calculate  the  mobility,  p.  We  obtain  the  following  expression  in  collision  time  approximation: 


_ e _ 

3ft2  kfiT 


~  ft)) 
Xkfo 


(10) 


The  mobility  calculated  by  including  the  scattering  due  to  ionized  impurities  and  polar 
optic  phonons  as  a  function  of  temperature  and  doping  density  is  shown  in  Figure  5a.  The 
change  in  the  Debye  screening  length,  complicated  variation  of  Fermi  level,  and  phonon 
scattering  give  raise  to  a  crossover  in  mobility  near  40  K,  as  seen  in  Figure  5a.  However,  it  is 
well  known  that  the  collision  time  approximation  always  overestimates  the  velocity  transition 
rate  that  determines  the  mobility. 

For  meaningful  comparison  with  experiment,  we  generalized  the  above  procedure  to  get 
full  solution  to  BTE.13  Thus,  calculated  mobility  with  our  hyperbolic  band  structure,  and  usual 
k*p  band  structure  is  shown  along  with  that  obtained  in  collision  time  approximation  and 
compared  with  experiments14  in  Figure  5b. 


One  important  feature  of  Figure  5b  is  that  the  hump  in  mobility  near  40  K  could  be 
explained  with  competing  impurity  and  phonon  scattering  rates.  Also  note  that  collision  time 
approximation  grossly  underestimates  the  mobility.  Although  k*p  values  are  closer  to 
experiment,  the  agreement  is  expected  to  be  poor  once  the  scattering  mechanisms  such  as 
acoustic  and  alloy  disorder  and  the  effect  of  temperature-dependent  band  structures  are  included. 
These  mechansism  lower  the  mobility  calculated  with  hyperbolic  band  structures  and  yield  an 
excellent  agreement  with  experiment 
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Figure  5a.  Mobility  with  hyperbolic  bands 
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Figure  5b.  Variation  of  Hall  mobility  with  temperature 
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CONCLUSIONS 


We  have  studied  the  effect  of  various  approximations  on  electron  transport  coefficients 
and  on  ways  to  extract  physical  parameters  from  experiments.  We  point  out  how  the  values 
interpreted  from  experiments  depend  crucially  on  various  approximations  such  as  effective  mass, 
MB  statistics,  and  collision  time.  The  results  are 

•  Approximating  the  Hall  factor  by  unity  over  a  wide  range  of  carrier  concentrations 
and  temperatures  is  accurate  for  most  applications.  An  error  of  about  30%  is 
expected  at  high  temperatures  and  low  carrier  concentration  (10 14  cm-3),  and  about 
20%  is  expected  at  low  T  and  high  carrier  concentration  (1018  cm-3). 

•  The  variation  of  band  gap  with  temperature  is  nonlinear  and  is  faster  at  low 
temperature.  The  value  at  22%  Hg  concentration  is  about  20  meV  smaller  than  values 
usually  quoted. 

•  The  mobility  calculated  from  full  solution  to  BTE  with  FD  statistics  can  explain  the 
hump  near  40  K  and  is  in  good  agreement  with  experiments.  Once  die  other 
scattering  mechansisms  such  as  acoustic  and  alloy  disorder  are  included,  the 
agreement  is  expected  to  be  much  better. 
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